Optical solitons in graded-index multimode fiber 
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Solitons are non-dispersing localized waves that occur in diverse physical settings. A variety of 
optical solitons have been observed, but versions that involve both spatial and temporal degrees of 
freedom are rare. Optical fibers designed to support multiple transverse modes offer opportunities 
to study wave propagation in a setting that is intermediate between single-mode fiber and free-space 
propagation. Here we report the observation of optical solitons and soliton self-frequency shifting in 
graded-index multimode fiber. These wave packets can be modeled as multi-component solitons, or 
as solitons of the Gross-Pitaevskii equation. Solitons in graded-index fibers should enable increased 
data rates in low-cost telecommunications systems, are pertinent to space-division multiplexing, and 
can offer a new route to mode-area scaling for high-power lasers and transmission. 
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Solitons are localized waves that arise from the inter- 
play of linear and nonlinear processes that individually 
would cause the wave to decay. They occur in numerous 
physical settings, including liquids [1, 2], optical fibers 
[3], plasmas [4], and condensed matter [5, 6]. Tempo- 
ral solitons that form in single-mode optical fiber [3, 7] 
are perhaps the quintessential example, and they have 
had major influence on telecommunications. The po- 
tential impact of soliton formation in multimode fiber 
was appreciated by early workers [8-10]. More recently, 
multimode fiber has been considered theoretically as an 
environment that could support spatiotemporal solitons 
[11, 12], or light bullets, which attract interest owing 
to their particle-like nature and potential for all-optical 
switching [13-15]. Perhaps surprisingly, no experiments 
have been reported despite the theoretical progress. 

Optical pulse propagation in a multimode fiber in- 
volves an intricate mix of spatiotemporal phenomena 
coupled fundamentally by nonlinearity and practically by 
waveguide imperfections. Modeling of the pulse propa- 
gation [8-10, 16, 17] is difficult, and multiple approaches 
have been taken. These include rewriting the coupled 
modes in terms of principal modes [18-20], variational 
solution of a nonlinear wave equation [11, 12, 21], and 
analysis of optical-wave thermalization and condensa- 
tion [22]. From a modal perspective, soliton formation 
requires nonlinear coupling between the modes to can- 
cel the effects of modal dispersion [8-10]. Such multi- 
component or vector solitons have been studied in several 
contexts [23-25]. There is only one prior report of spa- 
tiotemporal vector solitons [14], where the two compo- 
nents were different colors. Alternatively, nonlinear pulse 
pulse propagation in a multimode waveguide can be ana- 
lyzed with a three-dimensional Gross-Pitaevskii equation 
[26], which is widely used to model Bose-Einstein con- 
densates. The soliton solutions simultaneously balance 
nonlinearity with diffraction, dispersion, and a spatial 
harmonic potential. 

In addition to its intrinsic scientific interest, the for- 
mation of solitons in graded-index (GRIN) fiber will be 
relevant to applications. Owing to their low cost and 
ease of alignment, multimode fibers are widely used in 
high-speed local area networks [27]. The maximum data 
rate is limited by intersymbol interference that arises 
from modal dispersion. A new pulse propagation tech- 
nique that can retain the simplicity of multimode systems 
while avoiding modal dispersion should be beneficial to 
low-cost, high-speed systems. As systems approach the 
Shannon limit for information transmission, interest in 
exploiting multiple spatial channels [28-30], which could 
be transverse modes, has grown. Whether used to min- 
imize cross-talk in neighboring modes or to determine 
nonlinear limits to such approaches, soliton formation 
will be a factor in the design of such technologies. Fi- 
nally, there is a growing need for large-mode-area fibers 
for the generation and transmission of pulses with ever- 
higher peak powers [31]. Soliton transmission through 
large-core graded-index fibers could play a major role in 



these applications. 

In this Article, we describe theoretical and experi- 
mental observations of optical solitons and soliton self- 
frequency shifting in GRIN multimode fiber. Remark- 
ably, the stable solutions of the coupled-mode equations 
and the Gross-Pitaevskii equation are equivalent. Im- 
plications of the results for telecommunication, space- 
division multiplexing, and high-power laser and trans- 
mission systems are discussed. 




Figure 1: Soliton formation in GRIN fiber: Top: 
Schematic of the core and cladding of the fiber and the relative 
index as a function of radius. Middle: In linear propagation, 
the spatial modes separate in time, and individually increase 
in duration. Bottom: Propagation of the soliton. 



Results 

Coupled-mode analysis 

The most-complete theoretical model is based on cou- 
pled modes of the electromagnetic field [16, 17]. The 
complex electric field can be decomposed into a sum of 
spatial functions for the modes, F p (p, (/>, w), multiplied by 
the evolving envelopes, A p (z,lu): 

E{p, 0, w) = £ F p(P' h »>)e iMu) *M*> co). (1) 
p 

To simplify the problem, we consider only the radially- 
symmetric modes (written explicitly in the Supplemen- 
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tary Information), which are relevant to our experi- 
ments. The normalized modes, F p (p,uj), and correspond- 
ing propagation constants, f3 p (cu), can be written as [17] 
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where L p is a Laguerre polynomial, k = wno/c, A is 
the index difference between the center and the cladding 
of the fiber, R is the fiber core radius, and w = 
(2i? 2 /fc 2 A) 1 / 4 is the fundamental mode size. By ne- 
glecting higher-order dispersive and nonlinear effects, the 
equations can then be written as 
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where 5/3^ ) is the difference between the pth and 

the fundamental mode of the first (second) coefficient in 
the Taylor expansion of (3 p about uj . (3 2 corresponds 
to the material group- velocity dispersion, 7 = w «2/c is 
the nonlinear coefficient, and r/ plmn arc the nonlinear cou- 
pling coefficients (defined explicitly in the Supplementary 
Information) . The equations arc solved numerically for a 
variety of input fields launched into 100 m of a standard 
GRIN fiber (parameters in Methods), to correspond to 
experiments described below. The number of modes, and 
hence equations, required to account for the pulse prop- 
agation depends on the modes that are initially seeded. 
For the cases studied here the GRIN fiber is seeded by 
the output of a standard single-mode fiber and > 99.9% 
of the energy is accounted for with only the 3 lowest order 
symmetric modes (Figure 2(a)). 

In linear propagation, the modes separate in time ow- 
ing to their different group velocities and the pulse du- 
ration in each mode increases because of material dis- 
persion (middle panel of Figure 1). The beam waist os- 
cillates with period ^S== until the pulses separate tem- 
porally, after which the oscillation ceases and the spatial 
evolution is essentially that of the fundamental mode. At 
higher powers, nonlinearity balances the material disper- 
sion, mode-coupling counteracts the group- velocity mis- 
match, and a multi-component soliton forms. In other 
words, the duration and temporal location of light en- 
ergy in each mode are locked, and remain unchanged 
with propagation. While the existence of solitons in mul- 
timode systems was suggested based on analytical argu- 
ments [8-10], this is the first demonstration that multi- 
component solitons are even theoretically stable in a com- 
plete numerical model of the coupled modes. The result 
of launching a 300-fs, ^0.5-nJ pulse is shown in Figure 
2(b), e.g. The three modes overlap in time, with the 
group delay equal to an energy-weighted average of the 



modal delays. The spectra of the individual modes shift 
to have the same group velocity, with the higher-order 
modes blue-shifted (Figure 2(c)). The structured spec- 
tra result from radiation of energy from each mode as the 
soliton forms. Similar structure appears when solitons 
form in single-mode fiber (Supplementary Figure 2 in 
the Supplementary Information.) The space-time profile 
is nearly symmetric (Figure 3(a)), with a sech 2 tempo- 
ral intensity profile (Figure 3(b)) and a Gaussian spatial 
profile (Figure 3(c)). The mode-field diameter (MFD) av- 
eraged over the pulse (Figure 3(d)) oscillates around the 
fundamental mode size (2wo) with a period equal to that 
of low-intensity light. The full-width at half-maximum 
(FWHM) pulse duration averaged over the beam (Figure 
3(e)) quickly converges to a steady solution that is ~ 10 
times shorter than the output pulse duration would be 
due to group-velocity dispersion alone. 



Single-field analysis 

Alternatively, the system can be analyzed with a non- 
linear wave equation for the total field. Standard pro- 
cedures are used to reduce the Maxwell equations to 
a single wave equation. By employing the paraxial 
and slowly- varying envelope approximations, a Gross- 
Pitaevskii equation is obtained [11, 12]: 
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where A is the slowly-varying envelope at center fre- 
quency uj and t is time in a reference frame moving at 
the group velocity of the pulse. The paraxial approxima- 
tion does neglect the variation of group velocities among 
the modes. The approximation is justified in light of the 
solution to the coupled-field equations, which shows the 
modes aligning to a common group velocity. To reach 
manageable computational times, we reduce the system 
to the 2-dimensional case with y=0 as in Rcf. [12]. The 
solution quickly converges and simulations are stopped 
after 52 m of propagation (which requires 20 days of 
computation) . Remarkably, the soliton solution is nearly 
identical to that found with the coupled-mode equations 
(Figure 3), except for a slight decrease in the amplitude 
of the spatial oscillation (Figure 3(i)). 

The identification of terms in Equation 4 with specific 
physical processes provides insight, but exact analytic 
solutions are not known. An approximate analytic solu- 
tion is obtained by making the variational approximation 
with the trial function 
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Figure 2: Multi-component soliton: (a) The three spatial modes and the corresponding (b) temporal pulse profiles and (c) 
normalized spectra. The dotted lines correspond to the location of the pulse center after linear pulse propagation (without 
nonlinear attraction between modes). 




Figure 3: Simulation results: Numerical solution of the coupled-field equations (a-e) and the Gross-Pitaevskii equation (f-j). 
(a) Space-time intensity plot, (b) pulse averaged over space and (c) beam averaged over time (symbols) with a hyperbolic secant 
squared (solid line) and Gaussian (dashed line) fits, evolution of the (d) spatial averaged over the pulse and of the (e) FWHM 
pulse duration averaged over the beam from coupled-field equations and (f)-(j) the corresponding plots for the Gross-Pitaevskii 
equation. 



where E = J \ A\ 2 dxdydt is the energy, r is the pulse 
duration, wo is the beam width, 6 and a are chirp pa- 
rameters, and <f> is an arbitrary phase [11, 12, 32]. If the 

quantity e = 2 47r 2 |fe|fl IS mucn l ess than one, a stable 
fixed point to the equations of motion is given by 

w = (2i? 2 /fc 2 A) 1 /4 and E r = 2 ^ C7rw § . ( 6 ) 

uJon 2 

For typical fiber parameters e = 2 x 10~ 4 , so the approx- 
imation is excellent. The fixed point is thus the funda- 
mental mode of the GRIN fiber with a temporal profile 
that corresponds to the soliton with that beam size. 

Experiments 

The theoretical results suggest that the required pulse 
energy at 1550 nm can be reached with readily-designed 
Er-doped fiber lasers. Experiments were performed with 
a source that generates ~ 300-fs pulses (red lines in Fig- 
ure 4(a) and (b)) with energy up to ~ 3 nJ. The initial 



spatial profile is the fundamental mode of an ordinary 
single-mode fiber with MFD of 11.5 //m (Figure 4(c)). 
With these parameters, 100 m of GRIN fiber comprises 
~ 100 dispersion lengths (i.e., in linear propagation the 
pulse will broaden about 100 times), ~ 400 nonlinear 
lengths, and ~ 1,000,000 diffraction lengths. For en- 
ergies below 0.3 nJ the pulse disperses, and it is diffi- 
cult to measure the autocorrelation of the output pulse. 
Results of launching a 0.5-nJ pulse illustrate soliton for- 
mation (Figure 4). After propagation through the fiber, 
the pulse is compressed temporally (Figure 4(a)) and the 
spectrum becomes structured (Figure 4(b)). The MFD 
at the end of the GRIN fiber is 17.9 urn (Figure 4(d)). At 
higher energies, the output pulse duration decreases, and 
the spectrum broadens and red-shifts (Figure 5(a)), in a 
manner immediately reminiscent of soliton self-frequency 
shifting [33]. 

We can conveniently analyze the experimental results 
using the results of the variational approach. The depen- 
dence of pulse duration on pulse energy is predicted well 
by the fixed point (Figure 5(b)): the best fit is obtained 



with a MFD of 17.4 /im, which is close to the measured 
value of 17.9 ^m. For the spectral shift we apply the re- 
sults of standard soliton perturbation theory. The wave- 
length shift is inversely proportional to the fourth power 
of the pulse duration [33] : 
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Here z is distance and Tr is related to the slope of the 
Raman gain spectrum. The pulse duration is inversely 
proportional to the energy, so the wavelength shift should 
be proportional to the fourth power of the energy. Indeed 
this is the case (Figure 5(c)): the best fit is obtained with 
Tr = 2.6 fs, which is close to the accepted value of 3 fs 
for silica fiber. 
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Figure 4: Soliton formation with 0.5-nJ pulse energy: 

(a) Autocorrelation trace of the pulse and (b) spectrum before 
(red) and after (black) the GRIN fiber and beam profile (c) 
before and (d) after the GRIN fiber. 



Discussion 

The spatial and temporal degrees of freedom are 
coupled through the Kerr nonlinearity in a multimode 
waveguide. The picture that emerges from the theoret- 
ical and experimental results presented above is that a 
3-component soliton forms with nonlinear pulse propa- 
gation. That is, the process that is responsible for com- 
pensation of group-velocity dispersion in time also com- 
pensates for the separation of the different spatial modes 
(modal dispersion) . The spatial evolution of the solitons 
in GRIN fiber has a significant linear contribution, as 
expected with a waveguide. However, the system differs 
significantly from a single-mode waveguide because the 
spatial dimensions have the freedom to play a major role. 
Therefore, it is surprising that, on average (ignoring spa- 
tial oscillations), the resulting solutions are equivalent 
to single-mode solitons propagating in the fundamental 
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Figure 5: Experimental trends with pulse energy: (a) 

Output spectra for indicated pulses energies, (b) Measured 
pulse duration vs. energy (symbols) and fit to Equation 6. (c) 
Measured peak wavelength vs. pulse energy (symbols) and fit 
to Equation 7. 



mode of the fiber. This is noteworthy considering that 
the system includes multiple modes (Equation 3), and is 
modeled by a three-dimensional Gross-Pitaevskii equa- 
tion (Equation 4). 

Experimentally, measurements of the beam propaga- 
tion on a sub-millimeter length scale will be valuable to 
confirm the presence of multiple transverse modes, but 
we are unable to make such measurements with adequate 
accuracy. While the analyses presented above clearly 
account well for the experiments, we cannot absolutely 
rule out an alternative explanation. Although the models 
would become quite complex, Raman scattering could be 
incorporated directly, rather than treated as a perturba- 
tion. In such a model, the interplay of linear mode cou- 
pling (which distributes energy uniformly among modes) 
and Raman scattering could result in energy being "si- 
phoned" into the fundamental mode of the fiber [34-36] . 
The occurrence of this process would be noteworthy in 
its own right. We hope that the results presented here 
stimulate efforts to address this intriguing spatiotempo- 
ral problem. 

Regardless of the final theoretical interpretation, the 
observation of optical solitons and soliton self-frequency 
shifting in a GRIN waveguide has implications for a va- 
riety of scientific and technological issues. GRIN media 
offer a stable and convenient setting for the study of spa- 
tiotemporal wave propagation, with solutions that also 
pertain to Bose-Einstein condensates. Soliton transmis- 
sion in multimode fibers relaxes the stringent coupling 
requirements for single-mode systems, and thereby could 
reduce the cost of high-data-rate telecommunication sys- 
tems [27]. GRIN solitons will be pertinent to current 



■5 



efforts to reach the Shannon hmit [37] through space- 
division multiplexing [28-30]. In such systems, cross-talk 
can be minimized by utilizing modes that do not over- 
lap spatially as channels. Finally, the linear scaling of 
soliton energy and power with core diameter will benefit 
applications in high-power lasers and pulse transmission 
[31]. For example, commercially-available GRIN fiber 
with 1-mm core diameter will allow transmission of 100- 
kW solitons. 



Methods 

For the coupled equation simulations, A = 1550 nm, E 
= 0.5 nJ, n 2 = 3.2 x 10~ 20 m 2 /W, n = 1.444, /3 2 = -281 
fs 2 /cm, A = 0.0029, and R = 31.25 fim. An overlap- 
integral calculation with a 10- ^m Gaussian input reveals 
that > 99.9% of the energy is accounted for with 3 modes 
(92.22% in p = 0, 7.17% in p = 1, and 0.56% in p = 2). 
In linear propagation, the pulse in the p = 1 mode moves 
away from the pulse in the fundamental mode at a rate 
of 33 fs/m (99 fs/m for p = 2). 

We use a split-step Fourier-transform method to nu- 
merically solve Equation 4. With 1024 points in time 
and 512 in space, and a 100-^im step in the propagation 
direction, simulation of pulse propagation through 52 m 
of fiber requires 20 days of computation using four cores 
on an Intel i7 computer. The parameters used for the 
simulations are the same as those for the coupled system, 



with an effective length of 11.8 /xm for the y dimension. 

Experiments were performed at 1550-nm wavelength 
where the group-velocity dispersion of fused silica is 
anomalous. An Erbium-doped fiber oscillator (11 
MHz repetition rate) operating in the normal-dispersion 
regime produces ~ 300-fs pulses with <~ 1-nJ pulse en- 
ergy. The oscillator is followed by 80 m of normal- 
dispersion fiber stretcher, a single-mode amplifier, and a 
transmission grating compressor in a standard chirped- 
pulse amplifier configuration. The output is aligned into 
a SMF-28 fiber pigtailed collimator of 50-cm length, 
which allows for near ideal seeding conditions when it 
is spliced directly to the GRIN fiber (Thorlabs GIF625). 
After all losses due to the grating compressor and cou- 
pling into the collimator, this setup allows for up to <~ 3 
nJ at a dechirped duration of ~ 300 fs. 

The pulse is directly measured with a 2-photon in- 
tensity autocorrelator, the spectrum is measured with 
a grating spectrometer with 0.4-nm resolution, and the 
beam is measured with an InGaAs camera with a lllx 
beam expander to fill the detector array. The MFD is 
measured by averaging the widths of a Gaussian fit of 
the x and y cross-sections at the center of the beam. The 
pulse duration as a function of energy for the pulses after 
the GRIN fiber are measured by taking the full-width at 
half-maximum duration of an intensity autocorrelation 
and dividing by the correlation factor for a sech 2 inten- 
sity profile (1.543). 
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